Structural change in multipartite entanglement sharing: a random matrix approach 
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We study the typical entanglement properties of a system comprising two independent qubit envi- 
ronments interacting via a shuttling ancilla. The initial preparation of the environments is modelled 
using random-matrix techniques. The entanglement measure used in our study is then averaged 
over many histories of randomly prepared environmental states. Under a Heisenberg interaction 
model, the average entanglement between the ancilla and one of the environments remains con- 
stant, regardless of the preparation of the latter and the details of the interaction. We also show 
that, upon suitable kinematic and dynamical changes in the ancilla-environment subsystems, the 
entanglement-sharing structure undergoes abrupt modifications associated with a change in the mul- 
tipartite entanglement class of the overall system's state. These results are invariant with respect 
to the randomized initial state of the environments. 

PACS numbers: 03.65.Yz, 03.67.Mn,03.67.-a 
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Open-system dynamics involving environments com- 
prising only a finite number of elements have been proven 
a valuable arena for the study of interesting physical 
phenomena such as quantum chaos quantum ther- 
modynamics Q , entanglement and relaxation 0, 0] • In 
most of these studies, the use of random-matrix theory 
has proven technically very advantageous in modelling 
random collisions between parts of the overall system 
at hand. Random matrices have been helpful in deal- 
ing with many tasks of quantum information processing, 
including quantum data hiding, quantum state distinc- 
tion, superdense coding and noise estimation [5[. Very 
recently, random matrices have found extensive applica- 
tion in the characterization of Markovian decoherence ^ . 

In this paper we unveil a further interesting situation 
where the theory of random matrices finds fertile appli- 
cations: We study the typical amount of entanglement 
that can be set in a multipartite system comprising two 
environments of arbitrary (finite) size and a shuttling 
two- level ancilla that bridges their cross-talking (Fig. 1). 
Differently from Refs. [3|, \% [6|, |7|, random matrices are 
used in order to model the initial preparation of the en- 
vironments. These interact with the shuttling ancilla via 
a Hamiltonian model having pre-determined interaction 
strength. This allows us to investigate on the typical 
ancilla-environment as well as all-environment degree of 
entanglement simply by averaging the values correspond- 
ing to many hystories of random initial preparations. 
We point out the sensitivity of the entanglement-sharing 
structure on the kinematic and dynamical aspects of 
the interactions and its independence from the random 
preparation of the environments. In particular, we show 
that upon tuning of the coupling Hamiltonian regulating 
the ancilla-environment interactions, an abrupt transi- 
tion between two inequivalent classes of multipartite en- 
tanglement is achieved. The dimensions of the environ- 
ments and the number of interactions enter preponder- 



antly into the determination of the entanglement-sharing 
structure, as we quantitatively reveal. 

It is important to remark that our model does not al- 
low for phenomena of thermalization or homogeneization 
typical of effectively Markovian reservoirs such as those 
considered in Refs. In fact, the evolution here at 
hand is profoundly different from a forgetful open dynam- 
ics comprising environments which are weakly perturbed 
by the interactions with the ancilla. Memory effects are 
clearly seen in our analysis, which reveals how the state 
of each environment is deeply affected by its coupling to 
the shuttling ancilla. As we do not impose restrictions to 
the strength of each interaction, also the Born approxi- 
mation does not hold in our analysis. 

The remainder of this paper is organized as follows. 
In Sec. U we address the case of a single environment 
interacting with the ancilla A. Through numerical sim- 
ulations supported by a clear analytical study we show 
that when the ancilla interacts repeatedly with the same 
environmental qubit, the degree of entanglement is in- 
variant to both the random initial preparation and di- 
mension of the environment. Sec. [Ill extends our study 
to the case of two multi-qubit environments mutually 
connected by the bouncing ancilla. A repeatedly inter- 
acts (collides, as we often say in this paper) with two 
specific qubit s, each beloning to the respective environ- 
ment. In this case, the amount of entanglement that 
can be set is a delicate trade-off between both the di- 
mensions of the two environments and the order of the 
ancilla-environment interaction. This scenario allows us 
to address the main points of our study. First we show 
that, by tuning the form of the coupling Hamiltonian, 
the system undergoes an entanglement-sharing transi- 
tion between two non-equivalent classes of multipartite 
entanglement. As no time is explicitly involved, we re- 
fer to this effect as a "kinematic" transition. We then 
study how repeated sequential collisions are able to af- 
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feet the bipartite and tripartite entanglement in a way 
so as to mark a second, more dynamical entanglement- 
sharing transition. Sec. IIIII summarizes our findings and 
opens up perspectives for future work. Finally, an ap- 
pendix is devoted to the technical description of the way 
random unitaries are built in this work. 



I. SINGLE-ENVIRONMENT CASE: 
PROPAEDEUTIC RESULTS 

In this Section we study the single environment-ancilla 
system, which serves a useful case for the discussion of a 
few preparatory results that will then be reprised when 
the two-enviroment situation is studied. Let us consider 
a pure, initially separable state of an n-qubit environ- 
ment E = {ei, e2, .., e n }. As a reference-state for our 
analysis, we assume that all the qubits are initially in 
their respective ground state so that we have 

l^ = (g)|o> fe . (1) 

k=l 

We then introduce an ancillary two-level system, A, with 
which E interacts. The ancilla is assumed to be prepared 
in its ground state |0) A and it is meant to collide with a 
specific element of E, which we label e. The assumption 
on the preparation of A is simply matter of convenience 
as far as pure states are taken. Moreover, we will show 
that for dim (E) = 1 any state can be considered for the 
ancilla, including mixed states. The Hamiltonian model- 
ing such an interaction is taken to be of the anisotropic 
Heisenberg form 

3 

H A e = ^2 J ^ A ® °i,e (2) 




FIG. 1: (Color online). Sketch of the situation considered. 
An ancillary shuttle A, modeled as a two-level system, in- 
teracts with one of the elements of two spin environments, 
El,r, which are remotely located. The initial preparation of 
each environment state is random. Upon control over the de- 
tails of the interactions and the order of the collisions with 
A, effective manipulation of the entanglement-sharing struc- 
ture of the multipartite system comprising environments and 
shuttle is achieved. Environment Er is shadowed in order 
to indicate that we consider both the single-environment and 
two-environment cases. 



where Ji is an interaction strength and 0"; 5 A(e)' s are the 
Pauli spin operators of the A (e) qubit. We have used 
the labelling &\ = a x , &2 = cr y and &s = cr z . Our task 
is to investigate typical entanglement properties of the 
e — A system and we thus aim at removing any depen- 
dence of our analysis from the state of the environment 
E. In order to achieve this we proceed as follows: we 
prepare E in the state resulting from the application of a 
random n-partite quantum gate [9[ constructed using the 
random unitary matrix U r that we uniformly draw from 
the normalized Haar measure of the unitary group 
[l2[. In the Appendix we provide the technical details 
for the parameterization of such random unitaries, which 
are the building block of our numerical calulations. We 
then let the e — A system evolve as g = UhU r pU^U^ 

where U h = e~ iHAet and p = \tp) E (^| <g> \0) A (0| and de- 
termine the associated amount of entanglement. This is 
then statistically averaged over a large ensemble of ran- 
dom initial preparations so as to remove the dependence 
on the initial environmental state. 

We now anticipate a quantitative result achieved via 
our numerical calculations, which will then be justified 
in a fully analytical way: for an isotropic Heisenberg in- 
teraction defined by taking J x , y ,z = J in Eq. (j2j), the 
entanglement set between A and e only depends on the 
dimensionless interaction time Jt and [as far as the an- 
cilla is prepared in a pure state] is insensitive to the di- 
mension of dim(E'), and the preparation of A. One 
may think that this result is only due to the fact that A 
collides with a single environmental qubit, thus reducing 
our problem to a two-qubit interaction. However, this 
is definitely not the case. In fact, in general, the ran- 
dom evolution within the environment, encompassed by 
the random unitary matrices ?7 r 's, creates a multipar- 
tite entangled state. Although the ancilla A physically 
collides with e only, the state of the latter is crucially af- 
fected by the entanglement-sharing structure within 
which is highly non trivial. In Ref. [l3j], f° r instance, it 
was shown that the behavior of bipartite entanglement 
shared by any two elements of a multipartite register E 
is a decreasing function of dim(E). Intuition would then 
lead one to believe that the entanglement set between 
the interacting qubits would also depend on dim(£ T ) jl4j . 
Here we show that such a dependence is not in order. In 
order to understand this counter-intutive behavior, we 
take a two-step approach. First, by studying the case 
of dim(i^) = 1, we provide an intuition of the reasons 
behind the dependence of the typical amount of A — e 
entanglement on the sole interaction strength J. In this 
simple and yet useful case, a random unitary operation 
built according to the parameterization given in the Ap- 
pendix is represented by the 2x2 matrix 

Ur = e*( C ° S i e Z x Sin f"A (3) 
y — sin cpe %x cos <pe iyj J w 

where the angles a, ip ((/>) are picked up from the range 
[0, 2tt] ([0, 7r/2]) uniformly with respect to the normalized 
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Haar measure dU r — (87r 3 ) -1 sin (20) dfidipdxda [U, [HI, 
[l6| . In order to evaluate the entanglement set within the 
system, we use Wootters' concurrence [17]. For a general 
bipartite pure or mixed state described by the density 
matrix 0, concurrence is given by 

4 

C = max[0, ^/X~i~J2 ( 4 ) 

k=2 



where Ai > Xj (j = 2,3,4) are the eigenvalues of 
p (02 02) p* (02 (8> <r 2 ). When A is prepared in its 
ground state and upon evolution of the e — A system 
as g = UhUrpU^Ufr, we get the density matrix 



2iJt- 
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Jt 



c 2 W>) 

<i>+x)s(2Jt)S(24>) 
-W+x)c(2Jt)C(2ct>) 




-2iJt+i(*+x)5(2Jt)5(20) 
S 2 (</>)S 2 (2Ji) 
|S(4Jt)S 2 (0) 




-2ij*+i(v-+x)(7(2Jt) > 5(2^) 0\ 

-|S'(4Ji),S' 2 (</>) 

C 2 (2Jt)S 2 (4>) 

0/ 



(5) 



with C(x) = cos(x) and S(x) = s'm(x). It is then 
straightforward to see that the concurrence shared by 
e and A after the application of U r and Uh results in the 
elegant expression 

C eA = sin 2 ((/>)|sin(4Jt)|. (6) 

This shows that, given a specific preparation of the en- 
vironment, the only parameters governing the entangle- 
ment are (j) and Jt. The typical value of C e A is obtained 
by averaging the above expression over any possible uni- 
tary matrix U r uniformly drawn according to the proper 
Haar measure. Explicitly, we have to calculate 

~Y r 7 ^/^ 1 r^TT />27T p27V 

C eA = ^3 J sin(20) d(f> J dip J d\ J da C eA 
= i|sin(4Jr)|, 

(7) 

which may seem a special result arising from the pure- 
state preparation of the ancilla. However, this is defi- 
nitely not the case, as we now demonstrate. By starting 
with the mixed ancilla state 

pa=(* ,% o ), Poe[o,i] (8) 

and calculating the evolution of the e — A system, one gets 
a density matrix that is an easy generalization of Eq. (j5j) . 
By inspecting the eigenvalues of p (02 02) p* (02 ® 02), 
whose explicit form is too lenghty to be reported here, 
one finds no dependence on \-> i> or a i m full analogy 
with the case of Eq. ([7j). The explicit calculation of con- 
currence leads us to the expression 

C eA = + (2p - 1) cos(20)] sin(4Jt)|. (9) 

As J^ 2 cos(20) sin(2(j))d(j) = 0, the average concurrence 
(calculated using the appropriate Haar measure, as done 
before) turns out to be identical to Eq. (|6j). The study 



can be straightforwardly generalised to the case of qubit 
A being prepared in any coherent-superposition state, the 
only difficulty being a slightly more complicated expres- 
sion for the concurrence corresponding to any set prepa- 
ration of E. The message, however, is rather clear: re- 
gardless of the state into which the ancilla is prepared, 
the typical e — A entanglement is simply set by the 
rescaled interaction time Jt. This is strictly valid only in 
the statistical sense: if specific instances of preparation 
of both A and e are taken, such an independence does 
not hold anymore. However, contrary to a naive expec- 
tation, the typical entanglement does not vanish. As we 
see later, this result is the key to understand what occurs 
in the two-environment case. 

We now approach the second step of our proof by 
studying the invariance of the e — A entanglement with 
respect to dim(i?) when A is prepared in a pure state. For 
this task, we consider a simple extension of the previous 
case to a two-qubit environment E = {ei,e2}, initially 
prepared in a state described by 

1 3 

p E = -(l ei <g> l e2 +^/3fc,e 2 H ei ® 0fc,e 2 

k " 3 ( 10 ) 

+ y^ftfc,eid"fc,ei® ll e2 + ^ Xkl^k^ ® 07, e 2 ) 
k=l k,l=l 

with x the elements of the tensor accounting for the cor- 
relations between e\ and e 2 and (3 e . (j = 1, 2) the Bloch 
vector of qubit ej (j = 1,2) [I, |Q. This form holds 
for both entangled and separable two-qubit states and 
is thus a formal description of an arbitrary preparation 
of E. By taking e = e 2 (an arbitrary choice that does 
not affect the generality of our discussion) we follow the 
recipe for evolution described above. This time, before 
calculating the e — A concurrence, we have to trace over 
ei's degrees of freedom. Through a tedious but otherwise 
straightforward calculation, we see that although the tri- 
partite E — A density matrix depends on x an d (3 ei , the 
reduced density matrix of the e — A system only depends 
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on /3 e2 . By properly averaging over any possible prepara- 
tion of the environmental qubit we are led to a typical 
e — A entanglement that is identical to Eq. ((6]). These 
considerations can be extended to dim(E') = n, although 
the complexity of an analytical proof scales exponentially 
with n. However, our numerical calculations are in per- 
fect agreement with the analytic conclusions, supporting 
the general validity of the arguments used here and thus 
demonstrating the claimed invariance. 

The relevance of this result and its significance become 
evident when one considers that random-matrix theory 
is widely believed to provide a good effective descrip- 
tion of the state of a system in thermal equilibrium at 
a high temperature [l9j |. In such unfavorable conditions, 
the achievement of entanglement regardless of the initial 
state of the environment (i.e. its effective temperature) 
and its robustness against the complexity of E"s structure 
and the initial (pure) state of A is a remarkable feature. 
Our study reveals that through the use of a suitable in- 
teraction we can reliably create entanglement between a 
clean ancilla qubit and a mixed-state environment qubit. 
Furthermore, the invariance of the amount of generated 
entanglement with respect to the size of the environment 
strongly contrasts with the intuitive belief that C e A would 
have followed a monogamy constraints similar to those 
responsible for the bipartite entanglement in an n-qubit 
system [l3|, [20|. Besides its intrinsic interest, such an in- 
variance will be proven crucial for the understanding of 
the study conducted in Sec. El 



II. TWO-ENVIRONMENT CASE: 
ENTANGLEMENT-STRUCTURE TRANSITION 

We now turn our attention to the main scenario of our 
investigation, which consists of two independent and spa- 
tially separated environments, Er and El- This time, 
the ancilla A shuttles between them, colliding always 
with the same qubit of each environment. Again, we 
are interested in typical values of entanglement and, in 
order to cancel any dependence on the initial prepara- 
tion of the environments, we rely on a random-matrix 
approach. We first show that, in line with the results un- 
veiled in Sec.HJ one of the environments is entangled with 
A with a constant typical degree, dependent solely on 
Jt. This is accompanied by a few other results related to 
the details of the environment-ancilla interactions. Most 
importantly, however, we demonstrate that the genuine 
tripartite entanglement established among A and the en- 
vironmental qubits involved in the collisions is subjected 
to a drastic change in its structure, depending on the 
shape of the interaction Hamiltonian and the number of 
collisions. The fact that these results refer to typical (i.e. 
statistically averaged) quantities, make them even more 
intriguing. Not only multipartite entanglement can be 
generated: through suitable interaction engineering we 
can actually efficiently dictate its form and nature. 



A. Quantitative analysis of entanglement 

Let us first consider the case where dim(^) = 1 (we 
label its only element as {ei,}), while Er consists of 
qubits {ei, e2, .., e n }. If A collides only once, first with 
En and then El, we find that the degree of entanglement 
between the single-qubit environment El and the ancilla 
remains constant against the number of qubits belonging 
to Er, as show in Fig. [2] (a). We can easily understand 
this in light of the results discussed in the previous Sec- 
tion, where we have shown that the typical entanglement 
between A and a single-qubit environment does not de- 
pend on the input state of the ancilla. As the collision 
with Er, followed by the trace over its degrees of free- 
dom, simply prepares A in a mixed state, the invariance 
of C eL A against dim(.E#) is clear: quantitatively, C eL A 
is identical to Eq. (j6j). On the other hand, by calling 
cr the Er qubit with which A collides, we have that 
C eR A < C eL A, regardless of dim(ER), a behavior that 
can be rigorously explained as follows. For the sake of 
simplicity, we consider the simplest case where both the 
environments consist of a single qubit. By means of the 
decomposition in Eq. (|3]), we prepare a general (sepa- 
rable) initial state of cr and cl by specifying two sets 
of parameters {<fij,i/jj,Xj,ctj} (j = L,R). The ancilla 
is then let interact with cr and e/,, in this order. The 
degrees of freedom of cl are thus traced out so as to as- 
certain the concurrence between A and eR. This turns 
out to be 

C CrA = |sin 2 (^)cos(2Jt)sin(4Jt)|. (11) 

There is no dependence on the state of cl in this ex- 
pression, although the second collision has occurred and 
no average has yet been taken. However, this is quite 
clear in virtue of the fact that we are considering only a 
single interaction with each environment. After the col- 
lision with e^, there is no way for A to convey to eR 
information on the state of the left qubit. In fact, if we 
instead trace out cr, we find precisely Eq. 0. Look- 
ing for typical values, i.e. calculating averages over the 
tensor product of any possible unitary matrix by using 
the appropriate Haar measure introduced before, we get 
C GrA = |cos(2Jt)sin(4Jt)|/2 < C Gl a V Jt. Quantita- 
tively, for Jt = 1 the typical eR — A concurrence is 
c± 0.157, which perfectly agrees with the independent 
numerical estimate provided in Fig. [2] (a) (dashed line). 
The ordering relation involving C eR A and C eL A is pre- 
served as the number of qubits in Er grows, as shown 
in Fig. [2] (a), where the difference C eL A — Ce R A appear 
to increase with dim(£^^). This is a clear effect of the 
profound differences, stemming from the possibility to 
achieve multipartite entanglement, between the present 
situation and the single-environment case addressed in 
Sec. HI For dim(^) > 1, one cannot exclude that 
the random unitaries taken in order to prepare Er set 
some genuine multipartite entanglement among its ele- 
ments. Consequently, the entanglement set by the dy- 
namics studied here and shared by eR and A would be 
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(a) (b) 




dim(E R ) dim(Efl) 



FIG. 2: (Color online). Typical concurrence between the an- 
cilla A and the environment qubits €l,r for dim (.Ez,) = 1 and 
dim(ER) = 1,..,4. In panel (a), A interacts first with eR 
and then with cl, for a rescaled interaction time Jt — 1 (see 
inset). The blue dots (red squares) show the average concur- 
rence between A and eR (ez,). In agreement with the analytic 
formula in Eq. (0, C eL A = | sin4|/2 ~ 0.378. (b) Same as in 
panel (a), but with a reversed order of interactions. 

bound to follow many-body entanglement monogamy re- 
lations [20j, which lower its amount as Er grows in di- 
mension. This is also the reason behind the behavior 
shown in Fig. [2] (b), where by inverting the order of the 
ancilla-environment interactions we lose the invariance of 
the cl — A concurrence, an effect occurring for precisely 
the same reasons explained above. On the other hand, as 
before, the second interaction always sets the lowest de- 
gree of entanglement. It is worth remarking that Fig.[2]is 
the result of an independent numerical calculation that, 
although follows the formal recipe highlighted here, does 
not rely on the analytic result in Eq. (fTT|) . 

B. Entanglement-Structure Transitions 

The explicit introduction of multipartite entanglement 
in the process described until now directs us towards the 
main point of this study. Through the proper tuning of 
the dynamical and kinematic properties of the system, 
we shall demonstrate the ability to manipulate the shar- 
ing structure of the typical multipartite entanglement set 
among A, ez, and eR. For the remaining part of our study, 
we find it computationally more convenient to abandon 
concurrence for negativity [2l|, [22| . This is an entangle- 
ment measure devised from the Peres- Horodecki criterion 
on positivity of partial transposition (PPT). It is a nec- 
essary and sufficient condition for 2 x 2, 2x3 and ooxoo 
systems, and is a sufficient condition for an arbitrary sys- 
tem. For a bipartite case we define the negativity as [22j 

Mbi = max[0, -2\ neg ] (12) 

where X neg is the single negative eigenvalue arising after 
the partial transposition of the density matrix. The con- 
venience in using negativity stems for its straightforward 



generalization to the multipartite scenario: we simply 
have to consider all the possible bipartite splits in a given 
system and the geometric average of the negativity asso- 
ciated with any of them. This construction ensures that 
our measure remains an entanglement monotone [23[ . As 
we are mostly interested in at most three qubits, we con- 
centrate on tripartite negativity, which in our case reads 

Mri = [^A\e L e R K L \Ae R 3K R \Ae L }^ (13) 

where A^l^r-, for instance, indicates the bipartition of 
qubit A against the group of qubits and eR. 

In light of the analysis performed in Sec. Ill Al one can 
clearly conclude that, if tripartite entanglement is set 
within the e^ — A — eR system, this has to be W-class [24j . 
In fact, we find non-zero entanglement in any biparti- 
tion obtained by tracing out one of the components of 
the overall system. In Fig. [2j we have shown the en- 
tanglement between the ancilla and each environmental 
qubit involved in the collisions. Moreover, we have also 
checked that the — eR bipartition turns out to be in- 
separable in each of the situations addressed there. This 
is a property not shared by the GHZ-class [25] , where 
by tracing one qubit out of a tripartite register, one gets 
separable reduced states. In order to check that W-class 
entanglement is indeed in order here, we have utilized the 
witness [H, E3 

W=^1-|GHZ)(GHZ| (14) 

with |GHZ) = (|000) + |lll»/>/2 a three-qubit GHZ 
state and 3/4 the squared maximum overlap between 
a GHZ state and a W-class state. Eq. ([T4]) is a valid 
GHZ witness: a negative expectation value of VV suc- 
cessfully reveals GHZ-class entanglement. On the other 
hand, this construction guarantees that a positive expec- 
tation value is achieved for W-class states [26| . By gener- 
ating a sample of 1000 random environmental states and 
taking dim(£^^) up to 4 (so as to match the situation 
depicted in Fig. [2]), we have numerically verified that, as 
far as an isotropic Heisenberg model is used to model the 
collisions, (W) > for every single preparation of the en- 
vironments and regardless of dim(i^R). The calculation 
has been performed by maximizing the overlap between 
the sample states and |GHZ) over the tensor product of 
three local rotations, each acting on A, and eR re- 
spectively. Moreover, the tripartite negativity defined in 
Eq. (p~3|) turns out to be non-zero for any initial prepara- 
tion of our statistical sample, and so is the typical value 
obtained by averaging over it, which guarantees genuine 
tripartite entangled nature of the — A — eR states with 
strong evidence of their W-class nature. 

The task of this Section is to show that the typical 
entanglement-sharing structure of such tripartite state 
can be abruptly modified by biasing the coupling Hamil- 
tonian towards a specific spin non-preserving model. We 
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FIG. 3: (Color online). Bipartite and tripartite negativity re- 
sulting from a bilocal interaction model ruled by Eq. (|15[) , for 
single-qubit environments and |000) e ^ Ae ^ initial state. The 
full disappearance of entanglement at A = 1 occurs in virtue 
of this specific choice of initial state. At A = we have a point 
of entanglement-sharing change. Other three such points are 
visibile within the shown range of A's. 



thus recast Eq. ([2]) into (j = L, R) 



Hi 



b,j 



J[°1,J 



1 ^1, 



A QS> (Ti t 



(15) 



i=2 



where we have clearly taken J2,3 = A J (A G M) and J\ — 
J . The reasons behind this choice are best explained by 
means of the following analysis. For the sake of argument 
we focus on single-qubit environments and take qubits 
A, cl and 6r as prepared in their respective ground state. 
By considering an — A collision followed by an cr — A 
one, both ruled by Eq. (fT5|). we get 



-2i\Jt . 



, ehAeR = e —co8(Jt - XJt) |0) e Jcos(Jt - XJt) |00) 
- ism(Jt - XJt) |ll)] AeR -sin(Jt - XJt) \l) eL <8> 
[sm(XJt + Jt) |01) -icos(XJt + Jt) \10)] A e R - 

(16) 

Despite its innocence, the entanglement-sharing prop- 
erties of Eq. ([T6]) are quite interesting. In Fig. [3] we 
show the bipartite and tripartite negativity against the 
anisotropy parameter A, for Jt = 1 (arbitrary choice). 
The behavior of the bipartite entanglement is strongly 
dependent on the choice of A. In particular, at A = 0, 
which leaves us with H^^ = Jo\,a ® ^i,e j7 only two out 
of three possible bipartitions are inseparable: the e^ — eR 
subsystem remains separable regardless of the choice for 
Jt. This is easily checked by studying the partial trans- 
posed of the state resulting from \^) eLAeR upon trace 
over A. And yet, the tripartite negativity NtH is rather 
large at A = (cfr. Fig. [3l dashed line), demonstrating 
that we have a tripartite entangled state that, in evident 
contrast with the results of the previous Section, is out 
of the W-class. [28[. This qualitative feature is typical 
as it survives to a statistical average over random ini- 
tial preparations of the environments, their dimensions 
and the number of A — ej interactions (j = L,R). This 



latter feature is shown in Figs.[H where we plot the aver- 
age bipartite and tripartite negativity against A G [—5,5] 
and the number of collisions. Fig. [4] (c) shows that at 
A = there is no quantum correlations shared by the two 
environments, even in this statistically typical scenario. 

How can we understand this result and relate it to 
the claimed change in multipartite entanglement struc- 
ture? We refer to the studies conducted by Plesch and 
Buzek [29[ on the classification of multiparticle quan- 
tum correlations via entangled graphs. Following their 
lines, we represent a multi-qubit system with an inherent 
structure of shared entanglement as a connected graph. 
Each vertex embodies a qubit while a bond connecting 
two vertices represents bipartite entanglement shared by 
the components of the corresponding reduced state. For 
three qubits, four possible classes are identified, as shown 
in Fig. [5l In particular, Eq. (fTSjh as well as any of the 
states resulting from the application of random unitaries 
over the chosen reference state |000) e ^ Aei? , nave a sharing 
structure corresponding to the graph in Fig.[5](c), panel 
(a) being representative of W-class. As it is argued in 
Ref. [29], by following this graph-based classification of 
multipartite entanglement, a state showing tripartite cor- 
relations with the two-way residual entanglement struc- 
ture of panel (c) is GHZ-class. Moreover, a remark is 
due: such a two-way entanglement-sharing structure is 
possible, for pure states, only if there exist classical cor- 
relation between the separable qubits (e^ and cr in our 
case) . This property is readily confirmed in Eq. (fTBj) and 
any random element of our statistical sample. In fact, we 
find that 



+ Qe 



(17) 



where Q eL , R denotes the reduced density matrix of qubit 
£l,r after (random) preparation, unitary evolution and 
appropriate partial traces [30]. We can thus rightfully 
claim for the anticipated entanglement transition induced 
by the tuning of the effective order parameter embod- 
ied by A. The naturally generated W-class state arising 
from the use of an unbiased Heisenberg coupling is non- 
trivially changed into a GHZ-type entanglement charac- 
terized by a two-way sharing structure. The transition 
is typical in the sense explained so far: given any pure 
initial environmental state, we are able to create either 
type of tripartite entanglement structures simply by tun- 
ing the interaction properties. 

As anticipated, we also allow the ancilla to interact 
multiple times, in sequence, with each environment (each 
time with the very same cl,r qubit). Multiple interac- 
tions, however, do not appear to have any qualitative 
relevance on the generated entanglement structure, al- 
though there is an evident quantitative effect on the 
amount of entanglement shared both in bipartite and tri- 
partite sense. In Fig. [H (a) and (b) we plot the typical 
tripartite and bipartite entanglement of the A—eL system 
and it is easily verified that the same behavior is found for 
the A — eR system. For the whole considered range of val- 
ues for A the bipartite negativity is always non zero, sig- 
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FIG. 4: (Color online), (a) Tripartite negativity of the El — A — Er system against A G [—5, 5] and the number of collisions, 
for a dimensionless interaction strength Jt — 1 in an anisotropic Heisenberg coupling involving single-qubit environments, (b) 
Bipartite negativity for the A — Ej qubit-pair, regardless of j = L, R. (c) Bipartite negativity for the El — Er system. For 
A = there is no entanglement in this bipartition, while the ancilla-environment ones are inseparable. 



naling that this interaction model guarantees the genera- 
tion of entanglement between the ancilla and an environ- 
mental qubit even under the "minimum-control" condi- 
tions where one is unable to stop the ancilla-environment 
interactions after a desired number of collisions. More- 
over, as claimed above, at A = the e L — cr system is 
indeed separable, regardless of the number of considered 
collisions (see the M eL e R = line in Fig. [4] (c)). 




e L e R e L e R 

FIG. 5: (Color online). Entanglement graphs for the tripar- 
tite setting. Each sphere is a qubit, each bond represents a 
non-separable state, (a) W-class states having no separable 
bipartition. (b) Fully separable and GHZ-class states, with 
no residual bipartite entanglement, (c) Entangled-class with 
bipartition-dependent residual entanglement. If the state con- 
tains genuine tripartite entanglement, it can be enumerated 
in the GHZ class. Systems ez, and 6r share classical correla- 
tions, (d) Biseparable states encompassing a Bell pair [9[. 



It is very informative to study the effects that increas- 
ing the size of the environments has on the degree and 
type of shared entanglement for A = 0. We first take 
dim (Ex) > 1 and leave En with a single qubit. The 
scheme of interaction in the usual one, with being 
struck first. The features associated with the resulting 
entanglement sharing are shown in Figs. [6j In panel (a) 
we plot the tripartite negativity set within the state of the 
interacting qubits after tracing out the non-interacting 
elements from the El environment. We find that Aftri 
decreases with the dimension of Ex, which is a physi- 
cally meaningful result in light of monogamy relations 
that the system should adhere to [20]. A similar behav- 
ior is found for the bipartite ex — A entanglement, as 
shown in Fig [6] (b) , which can be easily interpreted in 
light of the discussion in Sec. Ill Al Finally, panel (c) 
shows the cr — A entanglement, which exhibits a con- 
stant trend that is consistent with our previous results 
on the isotropic Heisenberg Hamiltonian. We omit show- 
ing the 6r — e L entanglement as this is exactly null, in 
these conditions. The entanglement structure remains of 
the two-way GHZ class [29j although for dim(Ex) > 8, 
N eL A — 0. This implies the disappearance of the leftmost 
bond in Fig. [5] (c) but no further change, as genuine tri- 
partite entanglement is still present in the e L — A — eR 
system (see Fig. [6] (a)). 

The above analysis covers only partially the issue of 
increasing-sized environments, as it addresses an asym- 
metric situation. Another interesting effect arises if Er is 
let grow in size. In fact, the rate at which bipartite quan- 
tum correlations decrease in the case considered above 
suggests that, by increasing dim(E^), all residual en- 
tanglement would die off, eventually, leaving us with a 
sharing structure that reminds more of a typical GHZ 
state. In other words, we are looking for a configura- 
tion allowing an additional entanglement transition from 
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(a) (b) (c) 




FIG. 6: (Color online). Entanglement obtained for an increasingly-sized El environment (which comprises of up to eight 
qubits) and a single-qubit Er. We set A = in Eq. (|16p and the rescaled interaction time Jt = 1. (a) Tripartite negativity, 
(b) Bipartite eL — A entanglement against the dimension of El and the number of collisions, (c) Same as panel (b) but for 
the single-qubit environment Er. 



a graph of class (c) to one of class (b) in Fig. This is 
indeed possible. We consider dim^z,) = dim(^^) and 
investigate the behavior of both tripartite and bipartite 
negativity against the number of environmental qubits 
and the number of interactions with the ancilla. A typi- 
cal result, achieved by considering 18 pairs of collisions, 
is given in Fig. [7] (a) where it is evident that, by consider- 
ing larger environments, we are able to keep the tripartite 
entanglement non-zero using repeated collisions between 
A and the environments. Any two-qubit quantum cor- 
relation, on the other hand, disappears. The curves are 
the best-fits to the discrete numerical values obtained 
via our simulations. We have found that the functional 
form Ae~ B ^ dim ^ E ^~^ (with A, B e R) gives an excellent 
agreement with the numerical results. Fig. ([7]) (b) shows 

(a) 

Negativity dim 
0.8 I 

0.6 1\ 30 

0.4 ■ 1 \ 20 

0.2 I \ io 

Q.pr V ^. ■ . . u 

10 20 30 40 50 

FIG. 7: (Color online), (a) Negativity against the environ- 
ments' dimension dim(Ej) (j = L,R) for A = 0, 18 pairs 
of collisions and Jt = 1. The dashed line shows the tripar- 
tite negativity, while the solid lines are for the entanglement 
within each bipartite subsystem involving only one environ- 
ment. The eR — eL entanglement is obviously null. Each curve 
is rescaled with respect to its maximum value, (b): We plot 
the number of qubits required in the remote environments in 
order to make the tripartite and bipartite negativity smaller 
than 10 _ 6 when a fixed number of collisions is taken. The top- 
most curve is for Aftri while the other two curves, which are 
almost superimposed, show the cases associated with Af ej A- 



the number of qubits necessary in Ej (j = R, L) in order 
to get negativities smaller than 10 -6 (arbitrarily chosen) 
at a fixed number of collisions. It is clearly seen that 
the environmental dimension required to get entangle- 
ment lower than the chosen threshold is almost always 
much larger for the tripartite negativity than for the bi- 
partite ones, thus strengthening our conclusions: we have 
been able to dynamically drive the system state towards 
typical GHZ-class entanglement, therefore realizing yet 
another sharing-structure transition. 



III. CONCLUSIONS 

Using a random-matrix approach, we have shown that 
the typical entanglement properties of two remote qubit 
environments, indirectly communicating via the media- 
tion operated by a shuttling ancilla, can be effectively 
manipulated both quantitatively and qualitatively. Bi- 
partite as well as genuinely multipartite entanglement 
can be faithfully generated and shown to have intrigu- 
ing statistically averaged properties of resilience against 
arbitrary pure preparations of the environments, their 
dimension and the number of interactions that each has 
with the shuttle. Interestingly, by means of a simple 
tuning of the details of the interaction, we can achieve 
transitions in the entanglement-sharing structure from 
W to two-way GHZ-class [29]. We have studied the ef- 
fects that an increasing dimension of the environments 
have in such transitions. The use of random matrices 
has emerged, here, as a valuable tool for the agile han- 
dling of the numerical simulations required for the pur- 
poses of our study, enlarging the range of applicability 
of the technique to the context of multipartite entangle- 
ment generation under unfavorable conditions. This is 
a central problem in current quantum information and 
computation problems, where we hope that our results 
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will trigger further interest. We are currently investi- 
gating strategies, developed well along the lines of this 
paper, to study the interplay between quantum correla- 
tions and thermodynamical properties of a multipartite 
system. 
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APPENDIX 

Construction of Random Unitary Matrices 

Here we briefly outline the recipe to generate a random 
unitary matrix. The algorithm has been extensively used 
in recent works [ll|, [T2, l3ll - [33j . The parameterisation 
is based on the original work presented by Hurwitz in 
1897 (lOj . Any unitary matrix, U r , of dimension n, can 
be decomposed as 

U r = e ia E 1 E 2 ...E n _ 1 (A-l) 



where E{ is an n x n matrix. Matrices E^s are readily 
constructed using products of proper rotation matrices 
R^'^ (<t>i j i ipi j i Xi j) > eacn depending on the respective set 
of Euler's angles {4>ij,ipij,Xij} as follows 

E 1 =i? (1 ' 2) (012,^12,Xl2), 

E 2 = ^ (2 ' 3) (023,^23,O)^ 1 ' 3 )((/>i3,^13,Xl3), 
E S = ^ (3 ' 4) (034^34,O)^ 2 ' 4 )((/) 2 4,^24,O), 

X i? (1 ' 4) (0i4,^14,Xl4), (A-2) 

£ n _! = ^ n - 1 ' n )(0 n _ ln ,^n-m,O) X ... 

X # (1 ' n) (01n,^ln,Xln). 

The matrix elements are taken as 

Rff = cos 0, Rff = e ix sin 0, ( A-3) 

and zero otherwise. The angles are drawn from the 
ranges < faj < 7r/2, < ipij < 2tt, < Xij < 27r 
and < a < 2tt, uniformly with respect to the corre- 
sponding (and properly normalized) Haar measure 
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